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Abstract. Simple analytic arguments are used to understand the predominantly 
filamentary web in the large-scale distribution of galaxies. Numerical simulations 
of superclusters are performed to assess the feasibility of directly mapping the in- 
tracluster webbing, in particular with weak gravitational lensing using Wide Field 
Cameras. 

1 Patterns in random density fields 

In [[[] we have shown (see also [^) in this volume) that the large-scale structure 
in N-body simulations, which is believed to describe adequately the galaxy dis- 
tribution in redshift catalogues, can be explained directly by the geometry of 
the initial density field, amplified by gravitational instability. While rich clus- 
ters of galaxies form from high density enhancements of scales Rf ~ 8 h _1 Mpc, 
intracluster filamentary webbing reflects initial connections in the random den- 
sity field bridging these protocluster regions at lower density thresholds. 

The quantitative study of the properties of 3D random fields in a cosmolog- 
ical context goes back to the famous formula of Doroshkevich Q for the joint 
distribution function P(Ai, A 2 , A 3 ) of the ordered eigenvalues Ai < A 2 < A 3 of 
the initial deformation (shear) tensor ey. The latter is related to the density 
field, given in Fourier space by (random) amplitudes <5(k), by 

/r/ 3 k ~ - 
— kl k 3 S(k)W(kR f )e^. (1) 

The filter W(kRf) separates the smooth large-scale field which evolves in the 
mildly nonlinear single-stream regime from the highly nonlinear short-scale 
chaotic flow. 

Usually, the dynamical role of the shear in a local gravitational collapse 
is being focused on. Here we stress another aspect of the shear tensor: as a 
shape parameter describing the geometry of the initial density field. The mean 
overdensity profile around a field point ro with known values of the shear, 
(5(r)|Ai(r ), A2(r ), A 3 (r )), where 5 = p/p — 1, as illustrated in Fig. 1, shows 
that the shape depends on the sign of the shear eigenvalues (see Q for details of 
these and subsequent analytical calculations, also [||). In effect, the local shear 
signature defines the curvature of the density isocontours up to a distance of 
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Figure 1: Constrained primordial density field (5(r)|Ai, A2, A3) as a function 
of distance r in units of the filter scale Rf, in the three eigendirections (f-top, 
3-bottom). The left plot corresponds to the shear with all positive (+++) 
eigenvalues. The right plot represents the case of "filamentary" behaviour of 
the density in the neighbourhood of the (-++) sheared point. Dotted and 
dashed curves show the analytic short and long distance asymptotics [0 . 



several P/{] The increase in the density along one axis while falling off along 
the remaining two is symptomatic of filamentary bridges that connect the 
higher density regions where the shape of the density profile is more spherical. 
Fig. 2 shows that local properties of the shear tensor determine whether the 
region is part of the filamentary or protocluster structure. 

This correspondence between the global properties of the density field and 
local characteristics of the shear tensor enables analytical estimation of the 
gross properties of the Cosmic Web. For instance, the probability of the eigen- 
value signature given the density threshold 5, P(sign|J), shows which struc- 
tures dominates depends on the isodensity threshold chosen. Fig. 3 shows that 
at overdensities above a critical 5, 1.56a for Gaussian fields, one encounters 
predominantly spherical-like mass concentrations, while at the lower density 
contrast, < S < 1.56a, most of the initial density enhancements are in elon- 
gated filamentary bridges. Planar configurations ( h) are less likely for any 

positive overdensities S > 0. The related quantity P(<5|sign) gives us the den- 
sity distribution within different types of structure. While the average density 
of the filaments in the initial configuration is equal to 0.6a (Fig. 3), it is the 
1.5-2cr excursions which are precursors of the rare prominent filaments. In 
contrast, even rare planar, membrane-like, configurations are expected only at 
lower overdensities of 0.5-1 a. 

1 The information contained in the density curvature tensor 6 ij itself is much more local 
and less representative of the density behaviour at large distances from the constraint point. 



Figure 2: Example of the cosmological initial density field. Left panel: the 
part of space with overdensity 5 > la and A3 > A2 > Ai > 0. Note the roughly 
spherical shapes of the regions drawn. Right panel: the overdense regions with 
S > la and A3 > A2 > but arbitrary Ai. Inclusion of the points with Ai < 
creates filamentary connectors between round Ai > regions. 




Figure 3: Left panel: Probablity of the eigenvalue signature given the over- 
density threshold P(sign|i/), v = 5/a. Right panel: Density distribution given 
the type of shear tensor, P(V|sign). Mean densities for the given shear type 

are §d ~ 1.66er, 0.6a, — 0.6a for (+ + +),( — h +) and ( h) cases labeled by 

d = 3, 2, 1 with dispersion ad ~ 0.55ct almost equal for all configurations. 
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Figure 4: The nonlinear density evolution 5(t) of idealized spheres (top curve), 
filaments (middle) and membranes (bottom) of Lagrangian scale R / for which 
the smoothed linear rms amplitude of density fluctuations cr{Rf) is 0.65. Solid 
curves show 5(t) for the structures with initial (linear) density equal to 5d+2vd, 
d = f,2,3. The shaded bands delineate the range of 5(t) for initial densities 
[Sd + <Jd, $d + 2.5(7^]. Values of Sd and ad are given in Fig. 3. 



As a simple nonlinear model for evolution of clusters, filaments and mem- 
branes, let us consider the evolution of a shell which is spherical (d=3) or 
cylindrical (d—2) of radius R(t), or planar (d—1) of width 2R(t) which ex- 
pands with the Hubble flow in the other 3 — d directions (c.f. Fillmore and 
Goldreich 0|). The shell is characterized by its Lagrangian tophat scale Rf 
and the linear overdensity within the shell Sd- Evolution of the shell radius in 
different cosmologies with scale factor a(t) is given by Q 
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As we have shown, for cosmological initial conditions typical linear overden- 
sities 5d at the same scale Rf are quite distinct for structures of different 
geometry. Solving Eq.(j|) with Sd suggested by Fig. 3 demonstrates the range 
of present time densities one can expect for objects of different shape on the 
basis of this elementary theory. Results are presented in Fig. 4. For fila- 
ments we have p «3-f2p. Although filamentary structures have, generally, 
turned around, only rare filaments are approaching collapse by the present 
time. Membranes, which are still expanding, have p « 1-3 p. These values 
should be compared with the density of the collapsed spherical halos at the 
moment of virialization, « 180 p. N-body simulations confirm these simple 
estimates. 



2 flm and Ha are present-day density parameters of the matter and the cosmological con- 
stant. D(ti) is the growing mode at the onset of integration, R(ti) = a(ti) [1 — D(ti)S^/ d] Rf 



2 Observability of the Cosmic Web 



The filamentary pattern that joins galaxy clusters into a complex web is man- 
ifest in 3D galaxy catalogues. However the question of biasing between the 
galaxy distribution and the underlying mass density remains open. In re- 
cent work [j7j we evaluated the feasibility of mapping the intracluster struc- 
ture using observational techniques directly probing large-scale dark matter 
and gas distributions. We have performed gasdynamical simulations (using a 
SPH+TreeP 3 M code ||) of several rich supercluster regions and considered the 
weak gravitational lensing of high redshift galaxies by the low redshift struc- 
ture, the secondary variations in CMB temperature induced by the Sunyaev- 
Zeldovich effect and the X-ray emission by hot gas in clusters and along the 
filaments. The simulations and SZ maps derived from them are described in 
our companion paper 

Since the lensing depends upon the projected density E, not pressure as 
for the SZ effect or density-squared as for the X-rays, weak lensing would be 
a better probe for the filaments if the signal-to-noise is large enough. There 
are two sources of noise, one the intrinsic (random) ellipticitics of the source 
galaxies, and the other shear from ambient large scale structures along the 
line-of-sight confusing the ellipticity pattern created by the filamentary web. 

Results of our simulations for standard CDM, open CDM and a CDM 
model with a cosmological constant are presented in Fig. 5. They demonstrate 
that cluster-cluster bridges at redshifts z ~ 0.5 are sufficiently developed to 
produce lensing of background galaxies at the level of several percent which can 
be detected in the images of modern wide field cameras. One such detection 
of the bridge between two close clusters in the field of three at z = 0.4 is 
described in ||. 
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Figure 5: Colour maps of the convergence k = — 2E/£ cr it for our compact 
supercluster simulations at z = 0.2,0.5 and 1. E is the surface density (with 
the mean taken out), and the critical density E CT it depends upon the comoving 
distance to the lens and the source. We assumed the latter are at redshift z = 2 
for these plots. White areas correspond to n < —0.04 with grey contours 
marking the levels —0.08, —0.12, —0.16. Red dots denote high density regions 
k < —0.3. The level of convergence in the green areas is ~ —0.02. See Figs. 4,5 
of B for the corresponding dark matter distributions and SZ maps. 



